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Abstract 

The radiation loss of small-amplitude radially symmetric oscillons (long-living, spatially localized, time-dependent solu- 

tions) in two- and three-dimensional scalar field theories is computed analytically in the small-amplitude expansion. The 

amplitude of the radiation is beyond all orders in perturbation theory and it is determined using matched asymptotic 

series expansions and Borcl summation. The general results are illustrated on the case of the two- and three-dimensional 

sine-Gordon theory and a two-dimensional (j> 6 model. The analytic predictions are found to be in good agreement with 

H the results of numerical simulations of oscillons. 
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1. Introduction 

There is an increasing interest in long-living, spatially 
localized classical solutions in field theories - oscillons - 
exhibiting nearly periodic oscillations in time [H-[l(|. For 
the two and more spatial dimensional sine-Gordon the- 
ory, following the original naming in these objects are 



generally called pulsons [l7[-[2l|. Oscillons resemble the 
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"true" (time periodic and of finite energy) breathers of the 
one-dimensional sine-Gordon (SG) theory, but unlike true 
breathers they are continuously losing energy by radiating 
slowly. Just like a breather, an oscillon possesses a local- 
ized "core" , but it also has a "radiative" region outside of 
the core. Oscillons appear from rather generic initial data 
in the course of time evolution, in an impressive number 
of physically relevant theories including the bosonic sec- 
tor of the standard model [12]-[25||. Moreover they form 
in physical processes making them of considerable impor- 
tance [26l| — [33l] . 

A crucial physical characteristic of oscillons is the am- 
plitude of the outgoing wave determining their lifetime. In 
a previous work we have been able to perform an analytic 
calculation of the radiation amplitude of oscillons in scalar 
theories in the limit of small oscillon amplitudes in one spa- 
tial dimension (34[. In this limit one can perform an ex- 
pansion yielding breather-like configurations with spatially 
localized cores. The small-amplitude expansion yields an 
asymptotic series for the core, but misses a standing wave 
tail whose amplitude is exponentially small with respect to 
that of the core [3{| . Therefore in order to compute the ra- 
diative amplitude, more precisely its leading part, implies 
to go beyond all orders in perturbation theory. In the 
present paper building on the results of [II] we develop a 
method to compute the radiation amplitude in D = 2 and 
D = 3 spatial dimensions in the small-amplitude expan- 
sion. The main result of this paper can be summarized 
by the following simple formula determining the energy 



loss of a small-amplitude oscillon in D spatial dimensions 
(D < 4): 
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where Ad , Bry are given in Eq. (|4"5|) as functions of D and 
of the coefficients of the scalar potential. Although our 
analytic calculations are carried out in the limit of the os- 
cillon amplitude going to zero, the results arc valid for non 
infinitesimal values of the amplitude. We have found sat- 
isfactory agreement between the predicted energy loss Eq. 
(Q} and the "measured" one (by numerical simulations) 
for the SG theory in D = 2 and D = 3 and in a </> 6 -type 
model investigated in Ref. [34[ in D = 2. Our results also 
imply that no breathers depending on continuous param- 
eters can exist for D > 1 for arbitrary scalar potentials. 
Note that in D = 1 the SG theory is the only one with 
analytic potential admitting a breather family [36| . 

2. Small-amplitude expansion 

We consider spherically symmetric solutions of a real 
scalar theory in a D + 1 dimensional Minkowski space- 
time, with a self-interaction potential U(cj>). The equation 
of motion is 



dt 2 



d 2 4> 
dr 2 
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The mass of the field is chosen to be 1, and the deriva- 
tive of the potential U (</>) is expanded as a power series 
in (j>, where the Qk are constants. In the following un- 
less otherwise stated we shall consider potentials which 
are symmetric around their minima, i.e. g-ik = 0. 

The small-amplitude expansion of solutions of Eq. |(5J) 
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can be written as 
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where e -C 1 is the expansion parameter . Demanding the 
functions <f>k be bounded they all turn out to be periodic 
in time. For the class of symmetric potentials 4>2k = 0, 
k = 1,2,..., and the first two non-vanishing amplitudes 
are given as 
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P3 cos(wi) - 
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p\ cos(3wi) , 
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where A = —3^3/4, the frequency, to = y/1 — e 2 , and the 
functions p\ and p% depend only on the radial coordinate 
r. Introducing S and Z by 



Pi 



p 3 = fiZ. 



1 



1 



/' 



A ('24 



A 



8 



95 



((») 



the constants gk get eliminated from the equations deter- 
mining pi and P3, 



d 2 S D-l dS 



dp 2 p dp 
d 2 Z D-l dZ 
dp 2 p dp 
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where a rescaled radial coordinate has been introduced 
through p = sr. Eq. |(7J) admits regular solutions decay- 
ing exponentially in dimensions D < 4, implying similar 
properties for the regular solutions of Eq. (jHJ) . From now 
on we will only discuss the fundamental solution of Eq . 
(J7J without nodes. On solutions with nodes see Ref. [371 ] - 
The amplitudes (fofc+i, k > 2, are determined by linear 
inhomogeneous equations analogous to (J5J), and they all 
decay exponentially. The small amplitude series ^ does 
not converge in general - it is an asymptotic one. This 
is closely related to the fact that no true breathers ex- 
ist in generic theories because of the radiative degrees of 
freedom. Nevertheless the asymptotic series (0) yields an 
excellent approximation for the spatially localized "core" 
part of the oscillon (for details see Ref. (37l|). 

It turns out that the radiative amplitude of an oscil- 
lon is of the order of exp(— Bd/e) i.e. beyond reach of 
finite order perturbation theory in e. As found in Ref. [l2| 
any slowly radiating oscillon solution can be well approxi- 
mated by a suitable exactly time periodic, "quasibreather" 
(QB) having a well-localized core and an inevitable stand- 
ing wave tail, whose amplitude is in some sense minimal. 
In Section 4. the amplitude of this standing wave tail and 
that of the outgoing radiation will be directly related. To 
find this amplitude using the Segur-Kruskal method [35j], 
we have to analytically continue our functions to the com- 
plex plane and find the singularities closest to the real 
axis as a first step. Because ([5]) is linear in Z its solution 
has a singularity exactly where S has. Thus, we have to 



find the singularities of the master equation In the 
neighbourhood of the closest singularity to the real axis 
we compute the nonperturbative correction, which deter- 
mines the radiation field of the oscillon after matching it 
to the solution on the real axis. S has simple pole singular- 
ities on the imaginary axis at p = ±£Pd. Let us measure 
the distances from the upper singularity as 
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R. 
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Then the solution of {7J for small R can be written in the 
following series form 
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We note that for D > 1 the series (fTUf will also contain 
non-analytic terms, the lowest order being R 3 In R. The 
position of the pole can be determined both by Pade's ap- 
proximation and by integrating the master equation nu- 
merically along the imaginary axis. Table [T] gives the re- 
sults of these numerical approaches. For D = 4 Pjj = 0, 



D 


Pd 


1 


1.57080 


2 


1.09256 


3 


0.60218 



Table 1: The distance between the real axis and the pole of the 
fundamental solution of the master equation as a function of the 
number of spatial dimensions D. 

i.e. the pole is in the origin. It has been shown in Ref. 
[37( | that there are no regular and localized solutions of 
Eq. (|7|) for D > 4. This implies that in dimensions D > 4 
the oscillon core is not localized. As it will be shown the 
pole term Pd appears in the exponential in the radiation 
law, and from Table Q] it follows that oscillons in higher 
dimensions will radiate faster. 

The solution of the first inhomogenous equation (JS]) in 
the neighborhood of the singularity takes the form 



i2^2 1 
3 ~R? 



8V2(D 



15 P, 



D 



1) lni? z_ 2 



(11) 



The non-analytic terms represent dimensional corrections 
to the one-dimensional results. The free constant z_2 is 
due to the solution of the homogeneous part of the equa- 
tion for Z and is determined by the matching conditions, 
i.e. by using the fact that Z is decaying on the real axis. 

We have seen that the small-amplitude expansion yields 
time periodic and spatially localized (of finite energy) con- 
figurations which appear to be breathers. On the other 
hand any time periodic solution of Eq. can be expanded 
into a Fourier series as: 
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&2n = in the case of symmetric potentials presently con- 
sidered. For later use we compare the two different expan- 
sions, © and (H2): 
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Next we determine the behaviour of the Fourier modes, 
near the singularity. We define a new spatial coordi- 
nate y by R = ey, which is related to r as: 
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The "inner region" R <C f is not small in the y coordinates; 
if e — * then e\y\ = \R\ < 1 but \y\ — > oo. Then from 
Eqs. (fTT)|) and (jTTJ) the asymptotic behaviour for |y| — ► oo 
of S 1 and Z can be rewritten as 
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and from Eqs. (fT3"|) . we obtain 
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We shall impose Eqs. (|17[) . (fT5|) as boundary conditions 
to construct the inner solution (close to the singularity). 
We note that there is also a term proportional to elne/y 4 
in $3 arising from e 5 order terms in the small-amplitude 
expansion. 

To leading-order, i.e. setting e = in (JTTJ) an d (PT81) . we 
recover the one-dimensional result obtained in [34|]. The 
first dimensional correction comes from the elne terms. 
As we shall see, this will give a multiplicative correction 
factor to the energy loss rate. In the case of the SG model 
for D = 1 the breathers do not radiate. To determine 
the radiation loss in the SG models for D > 1 we have 
to calculate the next order corrections, arising from the e 
terms. 

3. Fourier mode expansion near the pole 

Let us now directly Fourier decompose our field equa- 
tion © using (|12p , without employing the small-amplitude 
expansion, 
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where the nonlinear terms on the right hand side have the 
form 



m,p,g— 1 



(20) 



Then in the inner region the mode equations (JT9J) keeping 
also 0(e) correction can be written as 



dy 2 
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By expanding equations (|2ip directly in powers of 1/y, 
including In y terms when necessary one can obtain Eqs. 
([T7|). (fig)) . Technically it is easier to obtain (TT7]) and (fT8"|) 
this way then by the small e expansion. 

Imposing Eqs. (JTTJ) , (fl8|) (and the corresponding ones 
for n > 3) as boundary conditions for Re y — > oo on the 
Fourier modes <!>„ defines a unique solution of the system 
(1211) . This corresponds to the behavior of the (real) so- 
lution for which all $„ decay exponentially when r —* oo 
on the real axis. Since in general no true breather exists, 
this decaying solution is singular at r = 0. Its extension 
to the complex plane is not real on the Re y = axis, 
the imaginary parts of the modes <&„ satisfy homogeneous 
equations to leading-order in l/y. Specifically, the solution 
for $3 behaves as 



Im$ 3 = v 3 exp(— iV8y) for Rcj/ = 0. 



(22) 



The constant v 3 determines the leading-order part of the 
radiation amplitude in the $3 mode which we aim to com- 
pute. A method to find the value of v 3 developed in Ref. 
[35j is to integrate Eqs. (f2Tj) along a constant Imy line 
numerically, starting from the values given by (fT7|). (fl"8|) 
for large \y\. An analytic method, using Borcl summation, 



has been proposed in Rcfs. 



and has been adapted 
for one-dimensional scalar oscillons in [3~j ] . 

Since to zeroth order in e Eq. (|2"Tj) is the same as for 
D = 1, the constant v 3 is independent of D for e — > 0, cor- 
responding to near threshold states with uj — ► 1. In order 
to find the corrections of order O(elne) and of 0(e) to 
the value of 1/3, it is sufficient to solve the mode equations 
linearized about the one-dimensional solution. Denoting 
the solution of the equations for D = 1 by &„ =1 + C( £2 )i 
and writing <!>„ = C E > .^ ,=1 + the equations linearized in 
e around D = 1 take the form 
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<3? n contains in general corrections of order 0(£m£) and 
of 0(e). The solution of the linearized equations to order 
C(£ln£) can be explicitly written as 



$ n = elneC— $^ =1 
dy 



(24) 



where C is an arbitrary constant independent of n. Eq. 
(I24p corresponds simply the zero mode when D = 1 in 
Eq. (|23p . The constant C is determined by the matching 
conditions (T7]) and (fTgj) : 



(25) 



In order to apply the Borel summation method described 
in detail in [34| , one expands the mode equations near the 
singularity, e.g. for <f> 



D=l. 
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and determines the growth of the coefficients Bk for large 
k as a first step. As found in Ref. [3]: 



B k ~ A' D=1 (-1) 



(2fc-2)! 
8 fc -!/2 



(27) 



The Borel sum of Eq. ([2^]) yields the imaginary part (in the 
case of symmetric potentials) of ^f =1 and by comparing 
it to Eq. (|22|) we have found the connection between the 
constant K and 1/3, yielding ^3° =1 = K D=1 tt/2. Expand- 
ing the mode equation near the singularity and calculating 
K from the behaviour of Bk for large k is a much easier 
and more precise way to calculate 1/3 than numerical inte- 
gration. Since the leading-order dimensional corrections, 
of order e In e are given by (|24|) , the dominant dimensional 
contribution to the value of K is given as 



K = K 



D=l 



1 - elneiVsC 
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where C is given in ([23|) . The relation ^3 = Kn/2 still 
remains true. In general, this expression already gives sat- 
isfactory results for the energy loss of oscillons, except for 
the SG theory, where K D=1 = 0. Then the correction of 
order 0(e) become essential. 

We now outline the computations of the dimensional 
corrections of 0(e) to the SG theory, in which case this 
correction is the leading one to the radiation amplitude. 
When K D=l 7^ there are other contributions to this or- 
der besides the dimensional correction, such as l/y correc- 
tions to Eq. ([22]) (see Eq. (23) of [li]) and the interaction 
of the outgoing wave with the core and these we do not 
compute here. It has already been pointed out that to each 
term of order 0(elne) there corresponds a term of order 
0(e), which is obtained by changing lne to \ny. Hence 
introducing $ n by 



$„ =elneC-^^ =1 +e fciny-^^ =1 
dy V d V 



the linearized equation (|23|) takes the form 
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The advantage of using is that there are no \ny terms 
in its 1 jy expansion. Expanding now $3 as 



— Pk 

— Z2k ' 
fc=l y 

from (|30p it follows that to leading-order 



(31) 



(2fc-2)(2fc-l)A_! + 8ft-(2fc-l)^B t = 0, (32) 

To find the large k behaviour of (3k, we consider the fol- 
lowing Ansatz: 



(33) 
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where i and M are the constants to be determined. Eq. 
yields L, 

r _ V2(D-1) rD=1 



8P D 



(34) 



while M is still arbitrary to this order in k. For the SG 
model K D=1 = 0, and by substituting the l/y expansion 
into Eq. (|30p we obtain a higher-order asymptotic formula 
for the behaviour of (3k for large values of fc, 



The value of M can then be obtained to a good precision by 
explicitly calculating the coefficients (3k in the expansion of 
(|30p to moderately high orders. Although we always work 
with a truncated SG potential and with a finite number 
of modes the value of K D=1 can be made very small by 
using <I>i — $9 and truncating the potential at 12 th order. 
For the constant K of the SG theory we obtain 



K = eM 



M = 0.6011 ; 



D-l 



P, 
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(36) 



We note that although a term proportional to z_2 (defined 
in (TTTj) ) also appears in the obtained values of (3k, its coef- 
ficient quickly tends to zero when increasing the order of 
truncation. Repeating the Borel summation argument for 
this case, it turns out that ^3 = Kn/2 still holds, giving 
the behaviour of Im $3 on the imaginary axis in the neigh- 
bourhood of the singularity. Naturally, the next step is to 
continue the imaginary correction back to the real axis. 

4. Continuation to the real axis 

The standing wave tail of a small amplitude time pe- 
riodic QB satisfies to leading order a homogeneous linear 
equation given by the left hand side of (fTT)]) . Since the 
size of a QB grows proportionally to 1/e, it is possible 
to consider an inner region, which is, however, outside 
of the domain around r = where the asymptotically 
decaying modes, $„, get large. In the previous section 
we constructed the function <5$3 by Borel-summing the 



4 



asymptotic series. We extend 5&3 to the inner region as- 
suming that it behaves as (|22[) close to the upper pole, 
and as Im$3 = — 1/3 exp(zv / 8?/) near the lower pole, where 
r = —iPd/s + V- The resulting function for large r is 



<5<I>3 = IC3 



(D-l)/2 



exp 



(37) 
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The general solution of the left hand side of (fT9|) can be 
written as a sum involving Bcssel functions J n and 3^, 
which have the asymptotic behaviour 



J v (x) 
Y v {x) 
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for x — * +00. The solution satisfying the asymptotics 
given by (|3"7) is 



*$ 3 = v^V^^r^-iCVSr) 



where the amplitude at large r is given by 
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This is just the solution singular at r = 0. The singularity 
is the consequence of the initial assumption of exponential 
decay for large r. The asymptotic decay induces an oscil- 
lation given by (|4"U|) in the intermediate core, and a sin- 
gularity at the center. In contrast, the QB solution has a 
regular center, but contains a minimal amplitude standing 
wave tail asymptotically. Considering the left hand side 
of (TH)|) as an equation describing perturbation around the 
asymptotically decaying solution, we just have to add a 
solution <5$3 determined by the amplitude (|4~Tj) with the 
opposite sign of (|4*0|) to cancel the oscillation and the sin- 
gularity in the core. This way one obtains the regular QB 
solution, whose minimal amplitude standing wave tail is 



given as 



= -^2^^^r>D/2-i(V8r)cos(3i) 

VSr-(D-l)^ cos(3i). 
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Adding the regular solution, where Y is replaced by J, 
would necessarily increase the asymptotic amplitude. 

Subtracting the regular solution with a phase shift in 
time, we cancel the incoming radiating component, and 
obtain the radiative tail of the oscillon, 
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The radiation law of the oscillon is easily obtained now, 
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(45) 

Without going into details in the case of asymmetric 
potentials, we notice that once the one-dimensional prob- 
lem is solved the first dimensional correction of order e In e 
is obtained in essentially the same manner, 

K = K D=1 (l-e In e 1 Vic) , A = h -g\ - \g z , (46) 
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The radiation law for small-amplitude oscillons is then 
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exp 



(47) 

Eqs. (05]) and (|4T|) can be interpreted as an evolution 
equation for oscillons of the corresponding theory. As 
the radiation power is exponentially small, the oscillons 
evolve through undistorted quasibreather states adiabati- 
cally. The oscillon energy was calculated in Ref. [13] and 



is E — e A ^ u ^ + e Ei as a function of e. Plugging 
this expression into the radiation law gives us the evolu- 
tion equation of e. 



5. Numerical simulations 

In this section we will provide numerical results regard- 
ing oscillon radiation. We have investigated a symmetric 
<fi 6 theory in D = 2, and the sine-Gordon theory in D = 2 
and D = 3. The results of numerical simulations in these 
theories confirm the theoretical predictions with satisfac- 
tory accuracy. The time evolution of oscillons was simu- 
lated using the fourth order method of line code developed 
in Refs. 40, 41 1. Initial data was obtained by applying the 
small-amplitude expansion method to e 3 order. A one pa- 
rameter fine tuning, by multiplying the amplitude with a 
factor close to one, was applied. For D = 2 this minimized 
the low frequency modulation of the oscillon state, while in 
D = 3 the tuning was used to suppress the single decaying 
mode. 

5.1. Oscillons in the symmetric </> 6 theory in D = 2 

In Ref. [3] we have found that oscillons of the sym- 
metric <f> 6 theory defined by the potential 



U{<P) = \4> 2 



I^ 4 + ^ 6 



(48) 



obeyed the theoretical radiation law in D = 1 to a sat- 
isfactory precision. In this theory the radiation of small- 
amplitude oscillons is maximal among symmetric 4> e the- 
ories, therefore we have a chance of observing the dimen- 
sional correction to the value of K. The e dependence of 
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Tabic 2: Radiation power W for the oscillons in the 6 theory. 
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Figure 1: Radiation law for small-amplitude oscillons in the 6 the- 
ory for D = 2. We plotted the theoretical radiation law with the 
first dimensional correction and without it. 



the energy loss rate W = dE/dt is listed in Tableland 
plotted on Fig. [TJ Analogously to the D = 1 case, we find 
anomalous points, corresponding to oscillons which radiate 
much slower than the theoretical prediction. This reveals 
the rich structure of finite e corrections and outlines the 
limitations of the calculation in the small-amplitude limit. 

Based on the theoretical radiation law (|4"5|) , in the D = 
2 case we use a semiempirical formula 



dE 

~dT 



-K 



,287.45 



■ exp 



6.18045 6 



(49) 



where an empirical correction constant 6 is included in 
the exponent, and the dimensional correction is included 
in the form of l[2"g]). with i\/8C = -1.592 in this case. 
The theoretical value of K D=1 has been obtained in [34j |. 
K D=1 — 0.579, while the theoretical value of 6 is ob- 
viously 6=1. We performed fits for the points with 
smaller e values than the last anomalous one. Fitting 



for both K D ~ Y and b for e < 0.25, we have obtained 
K D=1 = 0.05 ± 0.02 and b = 0.82 ± 0.02. We get the the- 
oretical pole term in the exponent quite precisely, but the 
amplitude K 1 has large error. However, setting 6=1 
and fitting only for K D=1 in the same domain, we have ob- 
tained K D=l = 0.59±0.05, which is close to the theoretical 
prediction. Interestingly, there are points outside the do- 
main of our fits which also obey the theoretical radiation 
law. It is evident from Fig. [JJ that the dimensional correc- 
tion significantly decreases the radiation power, but has 
little effect on the pole term, i.e. the curves have approxi- 
mately the same slope. The data points back the theoret- 
ical form of the dimensional correction, even though the 
precise e In e functional dependence is impossible to verify. 

5.2. Comments on the 4 theory in D = 2 

Oscillons in the cj) 4 theory show peculiar behavior in 
D = 2. In the e domain accessible to our lattice simula- 
tions they obey the semiempirical radiation law, however it 
is not consistent with the expected exp (— 2^/3P 2 /e) pole 
term, even though the outgoing radiation is dominantly in 
the $2 mode. We remark that in the </> 4 case the first di- 
mensional correction is so large in the accessible e domain 
that it invalidates the computed leading-order behavior. 

5.3. Oscillons in the D = 2 and D = 3 sine- Gordon theory 
In the D = 2 and D = 3 SG theory we expect oscillons 

to live exceptionally long, as the potential is symmetric 
and the dimension independent part of K is zero. Our 
simulation results are collected in Table [3l and plotted on 
Figs. H and H 
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Table 3: Radiation power in the SG theory. 

For the D = 2 case, we can again use the semiempirical 
formula (j49|) . where now K = eM. According to (|36|) . the 
theoretical value of M is M = 0.5502, and for the other 
parameter, obviously 6=1. First we fitted both M and 6, 
obtaining M = 1.79±0.16 and 6 = 1.08±0.01. Fitting only 
M while fixing 6 = 1 we got M = 0.76 ± 0.04. Similarly 
to the 4> e theory, the pole term can be measured quite 
precisely, but the amplitude has larger errors. 

Motivated by PS), in the D = 3 case we use the 
semiempirical formula 



dE 
~dt 



= —K' 



,190.81 



■ exp 



3.40644 6 



(50) 
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Figure 2: The radiation law for oscillons in the D = 2 SG theory. 
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Figure 3: The radiation law for oscillons in the D = 3 SG theory. The 
vertical line corresponds to the energy minimum at e = 0.342. Points 
to the left of the line correspond to stable oscillon states, while the 
others to states having one decay mode, requiring a one-parameter 
tuning. 



where K = eM. The theoretical values are M = 1.9964 
and b = 1. The fit for e < 0.342 yields M = 1.156 ± 
0.04 and b = 1.02047 ± 0.004. Again, the pole term is 
precise, but the amplitude has larger error. Setting 6=1 
and fitting only for M does not help in this case, it gives 
about half of the theoretical value, M = 1.0019 ± 0.01. 
Most probably this discrepancy is due the nonlinear effects 
related to the still too large values of e. Since the energy 
loss rate falls exponentially with e, our time evolution code 
cannot provide reliable results for smaller e values. 
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